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Q I Abstract 

We study the multiscale simplicial flat norm (MSFN) problem, which computes flat norm at various 
' scales of sets defined as oriented subcomplexes of finite simplicial complexes in arbitrary dimensions. 

We show that MSFN is NP-complete when homology is defined over integers. We cast MSFN as an 
I instance of integer linear optimization. Following recent results on related problems, the MSFN inte- 

i ger program can be solved in polynomial time by solving its linear programming relaxation, when the 

' simplicial complex satisfies a simple topological condition (absence of relative torsion). Our most sig- 

nificant contribution is the simplicial deformation theorem, which states that one may approximate a 
' general current with a simplicial current while bounding the expansion of its mass. We present explicit 

bounds on the quality of this approximation, which indicate that the simplicial current gets closer to the 
original current as we make the simplicial complex finer MSFN opens up the possibilities of using flat 
. norm to denoise or extract scale information of large data sets in arbitrary dimensions. On the other 

<^ ' hand, it allows one to employ the large body of algorithmic results on simplicial complexes to address 

. more general problems related to currents. 

o 

^ ■ 1 Introduction 

O ' Currents are standard objects studied in geometric measure theory, and are named so by analogy with 

electrical currents that have a kind of magnitude and direction at every point. Intuitively, one could think 
of currents as generalized surfaces with orientations and multiplicities. The mathematical machinery of 
currents has been used to tackle many fundamental questions in geometric analysis, such as the ones related 
r> ■ to area minimizing surfaces, isoperimetric problems, and soap-bubble conjectures [15]. We restrict our 

, attention to currents with integer multiplicities. 

The mass M{T) of a current d-dimensional current T can be thought of intuitively as the weighted d- 
dimensional volume of the generalized object represented by T. For instance, the mass of a 2-dimensional 
current can be taken as the area of the surface it represents. The flat norm of a d-dimensional current T is 
given by 

F(r) = mm{M(r - dS) + M{S) \ T-dS e^d, s e <^d+i}, (1) 

where is the set of d-dimensional currents with compact support and integer multiplicities. One also uses 
flat norm to measure the "distance" between two d-currents. More precisely, the flat norm distance between 
two d-currents T and P is given by 

F(r,P) = inf{M(Q) + M{R) \T - P = Q + OR, Q e <^d, R (2) 
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Morgan and Vixie |[T6]| showed that the total variation functional (L^TV) introduced by Chan and 
Esedoglu f3l computes the flat norm for boundaries T with integer multiplicity. Given this correspon- 
dence, and the use of scale in L^TV, Morgan and Vixie defined 1 16 1 the flat norm with scale A G [0, oo) of 
an oriented d-dimensional set T as 



where S varies over oriented {d + 1) -dimensional sets, and is the d-dimensional volume, used in place 
of mass. Figure [Uillustrates this definition. Flat norm of the ID current T is given by the sum of the length 
of the resulting oriented curve T — dS (shown separated from the input curve for clarity) and the area of the 
2D patch S shown in red. Large values of A, above the curvature of both humps in the curve T, preserve 
both humps. Values of A between the two curvatures eliminate the hump on the right. Even smaller values 
"smooth out" both humps as illustrated here, giving a more "flat" curve, as S can now be comprised of much 
bigger 2D patches. 



Figure 1: ID current T, and flat norm decomposition T — dS at appropriate scale A. The resulting current 
is shown slightly separated from the input current for clearer visuahzation. 

Figure [Uillustrates the utility of flat norm for deblurring or smoothing applications, e.g., in 3D terrain 
maps or 3D image denoising. But efficient methods for computing flat norm are known only for certain 
types of currents in two dimensions. For d = 1, Under the setting where T is a boundary, i.e., a loop, 
embedded in M? and the minimizing surface 5" € as well, the flat norm could be calculated efficiently, 
for instance, using graph cut methods llT3l - see the work of Goldfarb and Yin ifTTl and Vixie et al. [251, 
and references therein. Motivated by applications in image analysis, these approaches usually worked with 
a grid representation of the underlying space (M^). Pixels in the image readily provide such a representation. 

While it is computationally convenient that L^TV minimizers give us the scaled flat norm for the input 
images, this approach restricts us to currents that are boundaries of codimension 1. Correspondingly, the 
calculation of flat norm for 1-boundaries embedded in higher dimensional spaces, e.g., M^, or for input 
curves that are not necessarily boundaries has not received much attention so far. Similarly, flat norm 
calculations for higher dimensional input sets have also not been well-studied. Such situations often appear 
in practice - for instance, consider the case of an input set T that is a curve sitting on a manifold embedded in 
M^, with choices for S restricted to this manifold as well. Further, computational complexity of calculating 
flat norm in arbitrary dimensions has not been studied. But this is not a surprising observation, given the 
continuous, rather than combinatorial, setting in which flat norm computation has been posed so far. 

Simplicial complexes that triangulate the input space are often used as representations of manifolds. 
Such representations use triangular or tetrahedral meshes [9] as opposed to the uniform square or cubical 
grid meshes in and M^. Various simplicial complexes are often used to represent data (in any dimension) 
that captures interactions in a broad sense, e.g., the Vietoris-Rips complex to capture coverage of coordinate- 
free sensor networks [[5] [6]. It is natural to consider flat norm calculations in such settings of simplicial 
complexes for denoising or regularizing sets, or for other similar tasks. At the same time, requiring that 
the simplicial complex be embedded in high dimensional space modeled by regular square grids may be 
cumbersome, and computationally prohibitive in many cases. 



Fa(T) 



mm{Ya{T-dS) + XYd+i{S)} 



(3) 
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1.1 Our Contributions 



We define a simplicial flat norm (SFN) for an input set T given as a subcomplex of tlie finite oriented 
simplicial complex K triangulating the set, or underlying space Q.. More generally, T is the simplicial 
representation of a rectifiable current with integer multiplicities. The choices of the higher dimensional sets 
S are restricted to K as well. We extend this definition to the multiscale simplicial flat norm (MSFN) by 
including a scale parameter A. SFN is thus a special case of MSFN with the default value of A = 1. 

This discrete setting lets us address the worst case complexity of computing flat norm. Given its combi- 
natorial nature, one would expect the problem to be difficult in arbitrary dimensions. Indeed, we show the 
problem of computing MSFN is NP-complete by reducing the optimal bounding chain problem (OBCP), 
which was recently shown to be NP-complete [8 1, to a special case of the MSFN problem. We cast the prob- 
lem of finding the optimal S, and thus calculating MSFN, as an integer linear programming (ILP) problem. 
Given that the original problem is NP-complete, instances of this ILP could be hard to solve. Utilizing re- 
cent work [7 1 on the related optimal homologous chain problem (OHCP), we provide conditions on K under 
which this ILP problem can in fact be solved in polynomial time. In particular, MSFN can be computed in 
polynomial time when T is d-dimensional, and is (d + 1) -dimensional and orientable, for all d > 0. A 
similar result holds for the case when T is d-dimensional, and is (d + 1) -dimensional and embedded in 
M'^+i, for aU d > 0. 

Our most significant contribution is the simplicial deformation theorem (Theorem 15. Il l, which states that 
given an arbitrary d-current in \K\ (underlying space), we are assured of an approximating current in the 
(i-skeleton of K. This result is a substantial modification and generalization of the classical deformation 
theorem for currents on to square grids. Our deformation theorem explicitly specifies the dependence of 
the bounds of approximation on the regularity and size of the simplices in the simplicial complex. Hence it 
is immediate from the theorem that as we refine the simplicial complex K while preserving the bounds on 
simplicial regularity, the flat norm distance between an arbitrary d-current in [K| and its deformation onto 
the d-skeleton of K vanishes. 

1.2 Work on Related Problems 

The problem of computing MSFN is closely related to two other problems on chains defined on simplicial 
complexes - the optimal homologous chain problem (OHCP) and the optimal bounding chain problem 
(OBCP). Given a d-chain t of the simplicial complex K, the OHCP is to find a d-chain x that is homologous 
to t such that ||x||i is minimal. In the OBCP, we are given a d-chain t of K, and the goal is to find a 1)- 
chain s, of K whose boundary is t and |ls|li is minimal. The OBCP is closely related to the problem of 
finding an area-minimizing surface with a given boundary ||T51 . Computing the MSFN could be viewed, in 
a simple sense, as combining the objectives of the corresponding OHCP and OBCP instances, with the scale 
factor determining the relative importance of one objective over the other. 

When t is a cycle and the homology is defined over Z2, Chen and Freedman showed that OHCP is NP- 
hard 141 . Dey, Hirani, and Krishnamoorthy Q studied the original version of OHCP with homology defined 
over Z, and showed that the problem is in fact solvable in polynomial time when K satisfies certain condi- 
tions (when it has no relative torsion). Recently, Dunfield and Hirani [8| have shown that the OHCP with 
homology defined over Z is NP-complete. We will use their results to show that the problem of computing 
MSFN is NP-complete (see Section [2n ). These authors also showed that the OBCP with homology defined 
over Z is NP-complete as well. Their result builds on the previous work of Agol, Hass, and Thurston [1], 
who showed that the knot genus problem is NP-complete, and a slightly different version of the least area 
surface problem is NP-hard. 

The standard simplicial approximation theorem from algebraic topology describes how continuous maps 
are approximated by simplicial maps that satisfy the star condition ifTTl §14]. Our simplicial deformation 
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theorem applies to currents, which are more general objects than continuous maps. More importantly, we 
present explicit bounds on the expansion of mass of the current resulting from simplicial approximation. 
In his PhD thesis, Sullivan |22] considered deforming currents on to the boundary of convex sets in a cell 
complex, which are more general than the simplices we work with. But simplicial complexes admit efficient 
algorithms more naturally than cell complexes. We adopt a different approach for deformation than Sullivan. 
More significantly, our bounds on approximation are much smaller than the ones presented by Sullivan (see 
Section 15.21 ). Along with MSFN, our deformation theorem also establishes how the OHCP and OBCP 
could be used on general continuous inputs by taking simplicial approximations, thus expanding widely the 
applicability of this family of techniques to real-life data. 



2 Definition of Simplicial Flat Norm 

Consider a finite p-dimensional simplicial complex K triangulating the set 0, where the simplices are ori- 
ented, with p> d+1. The set T is defined as the integer multiple of an oriented d-dimensional subcomplex 
of K, representing a rectifiable d-current with integer multiplicity. Let m and n be the number of d- and 
(d + 1) -dimensional simplices in K, respectively. The set T is then represented by the d-chain X^™ ^ 
where ai are all d-simplices in K and are the corresponding weights. We will represent this chain by the 
vector of weights t G Z™. We use bold lower case letters to denote vectors, and the corresponding letter 
with subscript to denote components of the vector, e.g., x = [xj]. For t representing the set T with integer 
multiplicity of one, tj € {—1, 0, 1} with —1 indicating that the orientations of and T are opposite. But ti 
can take any integer value in general. Thus, t is the representation of T in the elementary d-chain basis of 
K. We consider {d + l)-chains in K modeling sets S representing rectifiable {d + l)-currents with integer 
multiplicities, and denote them similarly by SjTj in the elementary {d+ l)-chain basis of K consisting 

of the individual simplices Tj . We denote the chain modeling such a set S using the corresponding vector of 
weights s G Z". 

Relationships between the d- and {d+ 1) -chains of K are captured by its {d+ 1) -boundary matrix 
which is an m X n matrix with entries in {—1, 0, 1}. If the d-simplex is a face of the {d + l)-simplex Tj, 
then the (i, j) entry of [dd+i] is nonzero, otherwise it is zero. This nonzero value is +1 if the orientations of 
(Tj and Tj agree, and is —1 when their orientations are opposite. The d-chain representing the set T — dd+iS 
is then given as 

X = t - [dd+i]s. 

Notice that x G Z™. We define the simplicial flat norm (SFN) of T represented by the d-chain t in the 
{d + 1) -dimensional simplicial complex K as 

{m n I 

VVd(cT,)|:E,| + VVd+i(r,)|s,| I X = t - [dd+i]s, x G Z^ i . (4) 
i=l j=l J 

Since x and s are chains in a simplicial complex, the masses of the currents they represent (as given in 
Equation [B are indeed given by the weighted sums of the volumes of the corresponding simplices. The 
integer restrictions x G Z™ and s G Z" are important in this definition. It is unclear how one would interpret 
a fractional weight for a simplex in these chains, as we are studying integral currents. We generalize the 
definition of SFN to define a multiscale simplicial flat norm (MSFN) of T in the simplicial complex K by 
including a scale parameter A G [0, cxd). 

F^{T) = min \ Vd(ai)|xi| + A j Vd+i(r,-)|s,| | | x = t - [9d+i]s, x G Z^*" i . (5) 
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This definition is the simplicial version of the multiscale flat norm defined in Equation (|3]l. The default, or 
nonscale, SFN in Equation (U)) is a special case of MSFN with the default value of A = 1. Henceforth we 
study the more general MSFN. We assume the d- and (d + 1) -dimensional volumes of simplices to be any 
nonnegative values. For example, when cjj is a 1-simplex, i.e., edge, Vi(crj) could be taken as its Euclidean 
length. Similarly, Y2{Tj) for a triangle tj could be its area. For ease of notation, we denote Yd{o'i) by Wi 
and V(i+i(Tj ) by Vj, with the dimensions d and d + I evident from the context. 

Remark 2. 1. The minimum in the definition of MSFN (Equation |5]) indeed exists. The function 



is lower bounded by zero, as it is the sum of nonnegative entries (we have A > 0). Notice that Fg{T) = 
mins f^{T, S). Further, we only consider integral s defined on the finite simplicial complex K, and hence 
there are only a finite number of values for this function. Hence its minimum indeed exists, which defines 
the MSFN of t. On the other hand, the proof of existence of minimum in the original definition of flat norm 
for rectifiable currents employs the Hahn-Banach theorem lITOl pg. 367]. 

We illustrate the optimal decompositions to compute the MSFN for two different scales (A = 1 and 
A ^ 1) in Figure m Notice that the input set T, shown in blue, is not a closed loop here. It is a subcomplex 
of the simplicial complex triangulating i7. The underlying set Q need not be embedded in M? - it could be 
sitting in or any higher dimension. We do not show the orientations of individual simplices and chains 
so as not to clutter the figure. We could take each triangle to be oriented counterclockwise (CCW), with T 
oriented CCW as well, and each edge oriented arbitrarily. When scale A = 1, we get the default SFN of T, 
where the S chosen (shown in light pink) is such that the resulting optimal T — dS (indicated by the thin 
curve in dark green) is devoid of all the "kinks", but is similar to T in overall form. For A ^ 1, the second 
term in the definition (Equation[5ll contributes much less to the MSFN. As such, the optimal T — dS consists 
of a short chain of two edges (shown in light green), which closes the original T curve to form a loop. S 
in this case includes the triangles in the former choice of S, and all other triangles enclosed by the original 
curve T and the resulting T — dS. 



m 



n 



f\T,S) = Y, 



Wi\xi\ + X{2_^Vj\sj\) with X = t — [9(i+i]s 



(6) 




T - 95 for A «; 1 



Figure 2: MSFN illustrated for two different scales A = 1 and A ^ 1. See text for explanation. 



5 



2.1 Complexity of MSFN 



To study the complexity of computing the MSFN, we consider a decision version of the problem, termed 
decision-MSFN or DMSFN. The function f'^{T, S) used here is defined in Equation |6l with the modifica- 
tion that Wi and vj are assumed to be rational for purposes of analyses of complexity. 

Definition 2.2 (DMSFN). Given a p-dimensional finite simplicial complex K with p > d + 1, a set T 
defined as a d-subcomplex of K, a scale A € [0, oo), and a rational number /q > 0, does there exist a 
{d + 1) -dimensional subcomplex S of K such that f^{T, S) < /o? 

The related optimal homologous chain problem (OHCP) was recently shown to be NP-complete HI 
Theorem 1.4]. We reduce OHCP to a special case of DMSFN, thus showing that DMSFN is NP-complete 
as well. The default optimization version of MSFN consequently turns out to be NP-hard. 

Theorem 2.3. DMSFN is NP-complete, and MSFN is NP-hard. 

Proof. DMSFN lies in NP, as we can calculate f'^(T, S) in polynomial time when given a pair of d- and 
{d + l)-subcomplexes T and S, respectively, of the simplicial complex K. We just need to check that the 
chains t and s representing T and S satisfy x = t — [9rf+i]s, which can be done in time polynomial in the 
dimensions m, n of [S^+i]. Evaluating f^{T, S) can also be done in time polynomial in m and n. 

Conversely, we can reduce an equivalent instance of OHCP to a special case of DMSFN in a straightfor- 
ward manner. The input d-chain for the OHCP instance is t in the given simplicial complex K itself. Notice 
that the equations defining homologous chains in OHCP are given as c = t + [5rf_|_i]s, which are equivalent 
to the decomposition equations in DMSFN, as reversing the sign of s still gives a chain homologous to 
t. Consider the special case of DMSFN with A = and Wi = 1 for 1 < i < m. K solution to OHCP 
exists with the 1-norm of the homologous chain < /o if and only if this DMSFN instance has a solution 
/0(T, S) < /o. Since OHCP is NP-hard, we get that DMSFN is NP-hard as well. Therefore, DMSFN is 
NP-complete. 

A polynomial time algorithm for the optimization version of MSFN problem implies one for DMSFN. 
Given an instance of DMSFN, we apply this algorithm to the corresponding optimization version, and 
compare the flat norm value to /q. Hence the MSFN problem is NP-hard. □ 

We now consider attacking the MSFN problem using techniques from the area of discrete optimization. 
Even though the problem is NP-hard, this approach helps us to identify special cases in which we can 
compute the MSFN in polynomial time. 



3 MSFN and Integer Linear Programming 

The problem of finding the MSFN of the d-chain t (Equation |5]l can be cast formally as the following 
optimization problem. 

minimize ^™ ^ Wi\xi\ + KJ2]=i '^jl^jl) 

subject to X = t — [9^+1 ]s, C^) 
X G Z'", s G Z". 

The objective function is piecewise linear in the integer variables x and s. Using standard modeling tech- 
niques from linear optimization ||2] pg. 18], we can reformulate the problem as the following integer linear 
program (ILP). 

s.t. x+ - X- = t - [dd+i]{s+ - S-) (8) 
x+,x > 0, s+,s > 
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The objective function coefficients need to be nonnegative for this formulation to work - indeed, we have 
Wi,Vj, and A nonnegative. Integer linear programming is NP-complete fTF]. The linear programming 
relaxation of the ILP above is obtained by ignoring the integer restrictions on the variables. 



mm 



s.t. x+ - X- = t - [dd+i]{s+ - S-) 



x+,x > 0, s+,s > 



We are interested in instances of this linear program (LP) that have integer optimal solutions, which 
hence are optimal solutions for the original ILP (Equation |8]l as well. Totally unimodular matrices yield 
a prime class of linear programming problems with integral solutions. Recall that a matrix is totally uni- 
modular if all its subdeterminants equal —1, 0, or 1; in particular, each entry is —1, 0, or 1. The connection 
between total unimodularity and linear programming is specified by the following theorem. 

Theorem 3.1. || 24| Let A be an m x n totally unimodular matrix, and b € Z™. Then the polyhedron 
"P = {x E I Ax. = b, X > 0} has integral vertices. 

Notice that the feasible set of the MSFN LP (Equation |9l) has the form specified in the theorem above, 
with the variable vector (x+, x^, s+, s^) in place of x. The corresponding equality constraint matrix A has 
the form [/ — / B —B\ , where / is the identity matrix and B = [dd+i\- The input d-chain t is in place 
of the right-hand side vector b. In order to use Theorem 13. II for computing MSFN, we connect the total 
unimodularity of constraint matrix A and that of boundary matrix B. 

Lemma 3.2. If B = [c?(i+i] is totally unimodular, then so is the matrix A= \j —I B — -B]. 

Proof. Starting with B, we get the matrix A by appending columns of B scaled by —1 to its right, and 
appending columns with a single nonzero entry of ±1 to its left. Both these classes of operations preserve 
total unimodularity [ 18 , pg. 280]. □ 

Consequently, we get the following result on polynomial time computability of MSFN. 

Theorem 3.3. If the boundary matrix [d^+i] of the finite oriented simplicial complex K is totally unimodu- 
lar, then the multiscale simplicial flat norm of the set T specified as a d-chain t G ofK can be computed 
in polynomial time. 

Proof. The problem of computing the MSFN of T (Equation |5]) is cast as the optimization problem given 
in Equation This problem is reformulated as an instance of ILP (Equation [H). We get the MSFN LP 
(Equation |9ll by relaxing the integrality constraints of this ILP. As noted in Remark [ZTl the optimal cost of 
this LP is finite. The polyhedron of this LP has at least one vertex, given that all variables are nonnegative 
||2l Cor. 2.2]. By Lemma IJ!2l the constraint matrix of this LP is totally unimodular, as [dd+i] is so. Hence 
by Theorem 13. 11 all vertices of the feasible region of MSFN LP are integral, since t € Z™. 

An optimal solution (x+ , x~ , s+ , s~ ) of the MSFN LP can be found in polynomial time using an interior 
point method E Chap. 9]. If it happens to be a unique optimal solution, then it will be a vertex, and hence 
will be integral by Theorem 13.11 Hence it is an optimal solution to the ILP (Equation [H). 

If the optimal solution is not unique, then (x+, x^, s+, s^) may be nonintegral. But since the optimal 
cost is finite, there must exist a vertex in its polyhedron that has this minimum cost. Given a nonintegral 
optimal solution obtained by an interior point method, one can find such an integral optimal solution at a 
vertex in polynomial time [T2l. Hence the MSFN ILP can be solved in polynomial time in this case as 
well. □ 
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Remark 3.4. We point out that since the boundary matrix B = [dd+i] has entries only in {—1,0, 1}, the 
constraint matrix of the MSFN LP (Equation |9]l also has entries only in {—1,0, 1}. Hence the MSFN LP 
can be solved in strongly polynomial time [23 ], i.e., the time complexity is independent of the objective 
function and right-hand side coefficients, and depends only on the dimensions of the problem. 

Remark 3.5. Components of variables x+, x^, s+, in the MSFN ILP (Equation [8]) could assume values 
other than {—1,0, 1}, indicating integer multiplicities higher than 1 for the corresponding simplices in the 
optimal decomposition. The definition of MSFN (Equation [5ll does allow such larger multiplicities. At the 
same time, if one insists on using each {d+ l)-simplex at most once when calculating MSFN, and insists on 
similar restrictions on d-simplices in the optimal decomposition, we can modify the ILP such that Theorem 
E3]still holds. 

Denoting the entire variable vector by x = (x"*", x~, s"*", s~) G 22m+2n^ upper bound 

constraints x < 1, where 1 is the (2m + 2n)-vector of ones. These inequalities could be converted to the 
set of equations x + y = 1, where y is the (2m + 2n)-vector of slack variables that are nonnegative. These 
modifications give an ILP whose polyhedron is in the same form as described in Theorem 13.11 with the 
equations denoted as A'x.' = h' for the variable vector x' = (x, y). The new constraint matrix A' is related 
to the constraint matrix A of the original MSFN ILP given in Lemma [l!2l as 



A' 



A 6 
I I 



where / is the 2m + 2n identity matrix, and O is the m x (2m + 2n) zero matrix. Hence A' is obtained from 
A by first adding 2m + 2n rows with a single nonzero entry of +1, and then adding to the resulting matrix 
2n + 2m more columns with a single nonzero entry of +1. These operations preserve total unimodularity 
lITSl pg. 280], and hence the new constraint matrix A' is totally unimodular when [dd+i] is so. The new 
right-hand side vector b' S ^3m+2n consists of the input chain t and the vector of ones from the new upper 
bound constraints. 

Since the efficient computability of MSFN depends on the total unimodularity of the boundary matrix, 
we study the conditions under which total unimodularity of boundary matrices can be guaranteed. 



4 Simplicial Complexes and Relative Torsion 

Dey, Hirani, and Krishnamoorthy [T] have given a simple characterization of the simplicial complex whose 
boundary matrix is totally unimodular. In short, if the simplicial complex does not have relative torsion 
then its boundary matrix is totally unimodular. We state this and other related results here for the sake of 
completeness, and refer the reader to the original paper Id for details and proofs. The simplicial complex 
K in these results has dimension d + 1 or higher. Recall that a d-dimensional simplicial complex is pure if 
it consists of d-simplices and their faces, i.e., there are no lower dimensional simplices that are not part of 
some d-simplex in the complex. 

Theorem 4.1. [7, Theorem 5.2] The boundary matrix [dd+i] of a finite simplicial complex K is totally 
unimodular if and only if Hd{L, Lq) is torsion-free for all pure subcomplexes Lq,L of K, with Lq C L. 

These authors further describe situations in which the absence of relative torsion is guaranteed. The 
following two special cases describe simplicial complexes for which the boundary matrix is always totally 
unimodular. 

Theorem 4.2. [7, Theorem 4.1] The boundary matrix [da+i] of a finite simplicial complex triangulating a 
compact orientable {d + l)-dimensional manifold is totally unimodular. 
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Theorem 4.3. Q Theorem 5.7] The boundary matrix [dd+i] of a finite simplicial complex embedded in 
W^^^ is totally unimodular 

For simplicial complexes of dimension 2 or lower, the boundary matrix is totally unimodular when the 
complex does not have a Mobius subcomplex. 

Theorem 4.4. [7, Theorem 5.13] For d <1, the boundary matrix [dd+i] is totally unimodular if and only if 
the finite simplicial complex has no {d + l)-dimensional Mobius subcomplex. 

We illustrate the implications of the results above for the efficient computation of MSFN by considering 
certain sets. When the input set T is of dimension 1, and is described on an orientable 2-manifold to which 
the choices of 2-dimensional set S are also restricted, we can always compute its MSFN by solving the 
MSFN LP (Equation |9]) in polynomial time. A similar result holds when T is a set of dimension 2 described 
as a subcomplex of a 3-complex sitting in M^. For a 1-dimensional set T with choices of S restricted to a 2- 
complex K, we can always compute the MSFN of T efficiently as long as K does not have a 2-dimensional 
Mobius subcomplex. Notice that K itself need not be embedded in for this result to work - it could be 
sitting in some higher dimensional space. 



5 Simplicial Deformation Theorem 

The deformation theorem flO^ Sections 4.2.7-9] is one of the fundamental results of geometric measure 
theory, and more particularly of the theory of currents. It approximates an integral current by deforming 
it onto a cubical grid of appropriate mesh size. On the other hand, we have been studying currents or sets 
in the setting of simplicial complexes, rather than on square grids. Our proof is a substantial modification 
of the classical proof of the deformation theorem. We found the presentation of the latter proof by Krantz 
and Parks [14. Section 7.7] especially helpful. Our proof mimics their proof when possible. The gist of this 
theorem is the assertion that we may approximate a current with a simplicial current. 

Recall that Vrf((T) denotes the d-dimensional volume of a d-simplex a. The perimeter of a is the set 
of all its {d — 1) -dimensional faces, denoted as Per(cr) = {UjTj | Tj € fi, dim(rj) = d — I}. We will 
also refer to the {d — 1) -dimensional volume of Per(cr) as the perimeter of a, but denote it as P(ct) = 
St GPer(cr) Vd-i(rj). We let D(o-) be the diameter of a, which is the largest Euclidean distance between 
any two points in a. 

Theorem 5.1 (Simplicial Deformation Theorem). Let K be a p-dimensional simplicial complex embedded 
in M'^, with p = d + kfor k > 1 and q > p. Suppose that for every simplex a (z K 

D((j)P(ct) 
D(a) 

< K2 < OO, 

and 

D{a) < A 

hold, where is the largest ball inscribed in a, B^ is the ball with half the radius and same center as B^, 
and is the radius ofB^. Let T be a d-dimensional current in such that T C \K\, i.e., T is a subset of 
the underlying space of K. Suppose that T satisfies 

M{T) + M{dT) < OO. 
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Then there exists a simplicial d-current P supported in the d-skeleton ofK whose boundary dP is supported 
in the {d — \)-skeleton of K such that 

T-P = Q + dR, 



and the following controls on mass M hold: 



M(P) < (4^i^)^M(r) + A(4^x)*'+^M(9r), 
M(dP) < (4^i^)'=+iM(5T), 
M{R) < A(4^K)^M(r), and 
M(Q) < A(4^k)^(1 + 4dx) M{dT), 



(10) 

(11) 
(12) 
(13) 



where dx = Ki + K2. 

Remark 5.2. It is immediate that the flat norm distance between T and P can be made arbitrarily small 
by subdividing the simplicial complex to reduce A while preserving the regularity of the refinement as 
measured by ki and K2. 

Remark 5.3. Note that this theorem combines the unsealed and scaled versions of the original deformation 
theorem |[T4l Theorems 7.7.1 and 7.7.2] into one theorem through the explicit form of the constraints. In our 
proof of Theorem 15.11 we replace certain pieces of the original proof as presented by Krantz and Parks llT4] 
Pages 211-222] without reproducing all the other details of their proof. We found their exposition quite 
well-structured, making it easier to identify the modifications needed to get our theorem. 

Remark 5.4. The bound for M(P) in Theorem l5.1l is larger than the classical bound. We get this large bound 
because we generate P through retractions alone, and not using the usual Sobolev-type estimates 1141 Pages 
220-222]. And of course, the A in the coefficient of the extra term means that it becomes unimportant as 
the simphcial complex is subdivided. 

5.1 Proof of the Simplicial Deformation Theorem 

At the heart of the modification of the deformation theorem (from cubical grid to simplicial complex set- 
tings) is the recalculation of an integral over the current and its boundary. This integral appears in a bound 
on the Jacobian of the retraction, which measures the expansion in mass of the current resulting from the 
process of retracting it on to the simplices of the simplicial complex. To do this recalculation, we consider 
the retraction (j) one step at a time, building it through independent choices of centers to project from in 
every simplex and its every face. 

We first describe the general set up of retraction within a simplex. We then present certain bounds on the 
mass expansion resulting from the retraction in Lemmas 15.61 15.71 and 15.81 In particular, we obtain bounds 
on the expansion that are independent of the choice of points from which we project. These bounds are 
independent of the particular current that we retract on to the simplicial complex. But we employ these 
bounds to subsequently bound the overall expansion of mass of the current resulting from the retraction. 

5.1.1 Retracting from a center inside a simplex 

We describe the details of retraction for an ^-simplex a in the p-dimensional simplicial complex K. This set 
up is valid for any £, but in particular, we will use the bounds thus obtained ford<i<p when retracting a 
d-current onto K. We pick a center a G Int (cr), the interior of a, and project every x € Int (cr) \ {a} along 
the ray (x — a)/j|x — aj| to Fev{a). Denoting this map as i;^(x, a), we get 



</'(x,a) = ((/)^ o (/)5)(x,a) 



(14) 
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where (/>5(x, a) is a dilation of by the factor ||(/)(x, a) — a||/||x — a|| and (/)-n-{'^, a) is a nonorthogonal 
projection along (x — a)/||x — a|| onto Tx, the {£ — 1) -dimensional face of a containing (/)(x, a). We denote 
f = ||(/)(x, a) — a|| and r = ||x — a||. Let Ei be the ^-hyperplane that contains a and E^^i the {i — 1)- 
hyperplane that contains Tx. Denote the orthogonal projection of a onto Eg^i by b, and let /i = ||b — a||. 
For any point y = a + (b — a)7 with < 7 < 1, we get ^(y, a) = b. In particular, we consider the point of 
intersection of line connecting a and b with the {i — l)-hyperplane parallel to Tx that contains x. Naming 
this point y, we define ^ = ||y — a||. Let z € E"^ denote either normal to Tx at i?!)(x, a) (either of the two 
possibilities work). Let V2 = (x — a)/||x — a||, and let vi be the vector in span(z, V2) that is normal to V2 
and points into a. We illustrate this construction on a 3-simplex in Figure |3l where the cone of a with face 
r is shown in red and the other points and vectors are labeled. We also illustrate the corresponding slice 
spanned by vi and V2 in Figure |4] 




Figure 3: Illustration of the dilation and nonorthogonal projection involved in retraction for a 3-simplex. 



D0(<^(x))(vi) 




Figure 4: A 2-dimensional illustration of the dilation calculation. 
Choose an orthogonal basis {wi, w^_2} for {span(vi, V2)}"'". Note that span(wi, W£_2) G 
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Ei^i. Let w'be a unit vector in {span(wi, w^_2)}"'"n£'£_i parallel to 
V2} is an orthogonal basis for M^, and cf),^ is given by 



(x)-b. Then{vi,vi^i, ...,we^2, 



/i>7r(wi) 



aw 



Wj, i G {1, ...,£ — 2}, and 
0, 



(15) 



where a = f/h. Notice that the above set up works everywhere except when (/)(x) = b, in which case we 
obtain an orthogonal projection for (f).,^ (x) along b — a. Choosing coordinates for the tangent spaces of a 
and Tx to be {vi,-wi, vif£_2, V2} and {w', wi, ■Wi^2}, respectively, we get from Equation (1131 ) that 
1)0^ (x, a) is the {£ — 1) x £ matrix given as 



L>(/)^(x, a) 



a 
1 







1 



(16) 



5.1.2 Bounding the Integral of the Jacobian 

We now present a series of bounds on integrals of the dilation of d- volumes induced by the retraction. Since 
£ = d implies we are already in the d-skeleton and no retraction is needed, we can assume that £ > d. We 
start with a bound on the maximum dilation of d- volumes under the retraction cf). D(p will denote the tangent 
map or Jacobian map of 0. 

Definition 5.5. Let Jd'Pi'^, ^) be the maximum dilation of d-volumes induced by Dcl){-x., a) at x. 

We will use the definitions and results on D(/)7r(x, a) in ^-dimension given above. In particular, recall 
what D(cr), h, and h stand for. 

Lemma 5.6. For any center a and any point ain the £-simplex a with d<£<p = d + k, 

«(x,a)<(^-j -U. 

Proof. Following Equation (fT4)) . we seek bounds on Z)05(x, a) and ^^^^(x, a). Since D4>s{'^,sl) simply 
scales by r/r = h/h, the expansion of d- volume of any d-hyperplane by D(l)^{-x., a) is by a factor of {h/hY. 
On the other hand, bounding the dilation that L'07r(x, a) can cause in d-hyperplanes is a little more involved. 
We seek a bound on 

Vdet( (I)(/.^(x, a)^)^(Z)(/>^(x, aMT 
v'det(C/^C/) 

for all £ X d matrices U. Using the generalized Pythagorean theorem ||T4] Section 1.5], we get 

dei{U^U) = 5](det(C/A))2 

AeA 

where submatrix U\ consists of the d rows of U specified by the set of index maps A given as 

A € A = {/I/ : [1, d] [1, £], / is one to one and increasing}. 

A similar result holds for det((L>0^(x, a)C/)^(Z)0,r(x, a)C/)), with the functions / considered mapping 
[l,...,d]to [!,...,£ -I]. 
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Observe that multiplying by L'(/)7r(x, a) (Equation [T6l) just scales the first row of ?7 by a and removes 
the last row. Thus adet{Ux) > det((L'(/)7r(x, a)C/)A), which implies that a is a bound on the ratio in 
Equation ([TT] ). Thus we have that 



Jd0(x,a) 



r\'^ ||0(x,a) 




|b-a|| 

holds for all x and a in a, where D(cj) is the diameter of the ^-simplex a. □ 

Next we describe a bound on the integral of Jd4>{'^, a) over the entire ^-simplex, for a fixed center a. 
We will find that this bound is independent of the position of a. Recall that Per{a) and P{a) denote the 
perimeter of ^-simplex a and the {i — 1) -dimensional volume of the perimeter, respectively, and lnt{a) its 
interior 

Lemma 5.7. For any fixed center a in the i-simplex a with d<i<p = d + k, 

Jrf,/.(x,a)d/:^(x) < D(a)P(a). 



Int((T) 

Proof. Consider the {£ — 1) -dimensional faces tj of a, with Per(o") = {UjTj \ tj G a, dim(Tj) = £ — 1}. 
Let aj denote the ^-simplex generated by a and tj. Then 



Jd0(x, a) d/:^(x) = Y, [ Jd4>{x, a) dC\^) 

Int((T) J Jlnt(CTj) 



Let Tj{h) denote the {£ — l)-simplex formed by the intersection of aj and the (I — l)-hyperplane parallel 
to Tj at a distance h from a. Thus, Tj{h) is r itself. We observe that our bound on Jd4'{'^, a) is constant in 
Tj{h) for any h. The [i — 1) -dimensional volume of Tj{h) is given by 



Using the bound on Jd4>{^, a) from Lemma ls!6l and noting that D(cjj) < D(cr) V j, we get 



[ J,0(x,a)d£^(x)< / (i) " V,_i(r,) ^dh 
Jintfa,) Jo \hj yn h 



d 



'lnt(a-j) 

Summing this quantity over all Tj G Per(o") and replacing £ — d>l with 1 gives the overall bound. □ 

We now bound the integral of Jd4>{'^, a) over centers a with a fixed x that we are retracting onto Per((7). 
Examination of the corresponding proof for the original deformation theorem llT4l Section 7.7] shows that 
symmetry of the cubical mesh plays a very special role, which cannot be duplicated in the case of simplicial 
complex. In particular, we must avoid integrating over a close to the perimeter of a. Hence we integrate 
over as big a region as we can while still avoiding a neighborhood of the perimeter As in the statement of 
the main Theorem 15.11 let S^r be the largest ball inscribed in a, Bfj be the ball with half the radius and same 
center as B^, and be the radius of B^j. 

Lemma 5.8. For any point x in the l-simplex a with d<£<p = d + k, 

B(a) 



[ Jrf,/.(x,a)d/:^(a) < D(a)P(a)+V,(S,)^^ 
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Proof. Similar to the subsimplices of a considered in the Proof of Lemma 15.71 let aj now denote the £- 
simplex formed by x and Tj G Per(cr). In order to derive an upper bound, we integrate instead over regions 
that are by construction bigger than these subsimplices of a. Denoting the simplex aj as Region 1, we define 
Regions 2 and 3 as follows. We refer the reader to Figure 2]for an illustration of this construction. Let be 
the reflection of aj through x, and similarly, let rj be the reflection through x of Tj. We define a'- as Region 
2. Notice that unlike Region 1, Region 2 need not be contained fully in a. As defined in Section [5. 1.1 [ let 
z be the unit vector normal to the {i — l)-hyperplane containing Tj pointing into a. We define Region 3 as 
the ^-dimensional set rj + [0, D((t)]z, as illustrated in Figure |5] 

Note that the union of all Region 2's and Region 3's cover a. By an argument almost identical to that 
above, we have the following upper bound on the integrand in question. 



h' 



(x,a) 



< 



h 



h'YT){a) 



h 



Integrating the second of these two terms over Region 2 and summing the integral over all such Regions 




D(S) 



Figure 5: Illustration for integration of Jacobian bound over centers a instead of x's. For the case of the 
triangle shown, there will be 3 sets of 3 regions. In general there will be 3 regions for every face of the 
simplex. 

2 for all faces Tj, we get the upper bound of D((t) P(ct). Here we use the same arguments as the ones 
employed in Lemma [5^ Region 2 alone is not guaranteed to cover as some of Ba may occupy parts of 
Region 3. Since a G Scr, we have h > r^, and h' < h when a GRegion 3, so that 

(h'^Bja) ^B{a) 
\hj h ~ ' 
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Combining the above estimates while integrating over all such Regions 2 and 3 gives us the bound specified 
in the Lemma. □ 



5.1.3 Bounding the pushforwards of the current 

We consider the d-current T, and employ the bounds on the Jacobian of retraction described above to the 
pushforwards of T and its boundary dT on to the simplicial complex K. Our treatment of the pushforwards 
essentially follows the corresponding results of Krantz and Parks for the case of square grid llT4] Pages 
218-219]. We denote by ||T|| the total variation measure of the current T. 

Lemma 5.9. Suppose K is a p-dimensional simplicial complex with p = d + k for k > 1. Consider 
the stepwise retraction of the d-current T C K (the (d — l)-current dT C K) onto the d-skeleton of K 
(respectively, the [d — l)-skeleton of K). Each step of the retraction on to the perimeter of an i-simplex a 
for d < I <p (respectively, d < i < p) increases the mass ofT or dT by at most a factor of 

4.dK = 4(ki + K2) = 4max ; ' \^ + ^ ^ 



Proof. Using Fubini's theorem [ 14, Page 26] and applying the bound in Lemma |5^ we get 

[ /" Jrf(/)(x,a)d||r||(x)d/:^(a) = / / Jd0(x,a)d/:^(a)d||r||(x) < ^,M(ru), 

where -9^ = D((t)P((t) + N ^(Bff){T>{(j) Ir^j) and Tl^ is the portion of the current T restricted to the 
simplex a. Consider the subset of defined as 

/^T = |a G I ^ «x,a) d||r||(x) > l^^^gMj . 

Then \^{Ht) < (1/3) \ i{Ba-)- Similarly we define Hqt for the pushforward of dT and get Y^{HffT) < 
(1/3) V£(ficr). Then the set B(j\{Ht^Hqx} defines a subset of with positive measure, with the centers 
a in this subset satisfying Jd0(x, a) d||r||(x) < Ad^M{T\^) / \ ^{B,) and Jd0(x, a) d||5r||(x) < 
4^CT '^{T\a) I 'Vi{Ba)- Hence we can choose centers to retract from in each simplex a such that the expan- 
sion of mass of the current restricted to that simplex is bounded by AdafVeiBa-)- The bound specified in 
the Lemma follows when we consider retracting the entire current over multiple simplices in K, and set 
•&X = max^ex as the generic upper bound that holds for all simplices in K. □ 



Bound on complete sequence of retractions. We can apply the bound specified in Lemma 15.91 over 
multiple levels i. Pushing T onto the d-skeleton of p-complex K multiplies the mass of T by a factor of at 
most (4^^:)'^. Likewise, pushing dT on to the (d — l)-skeleton multiplies the mass of dT by a factor of at 
most (4^^-)'=+!. 



5.1.4 Bounding the distance between the current and its simplicial approximation 

In the final step, we construct the simplicial current P approximating the original current T, and bound the 
flat norm distance between the two. Since we are now considering retraction maps over many simplices 
simultaneously, we refer to all such maps as (p in general, suppressing the particular x and a. Hence we map 
T forward in k applications of 0, picking centers (see Lemma 15^ to project from in each step and in each 
simplex. We denote the composition of all these steps as tp^. We pick each of these centers such that the 
retractions map dT with bounded amplification of mass as well (see Lemma ls!9l ). 
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Define the homotopy (7(7, x) = 7X + (1 — 7)^^(x) for 7 € [0, 1]. Then the homotopy formula |[T4l 
Section 7.4.3] gives 

T = ^1{T) + dg#{[0, 1] X T) + 5#([0, 1] x dT). 
We define R = g#{[0, 1] x T) and Qi = g#{[0, 1] x (9r). Then we get 

T-i^l{T) = dR + Qi. (18) 

Finally, we map ip^idT) forward to the (d— 1) -skeleton of simplicial complex K with one more application 
of to get tp^idT) = (/'#(V'#(9T)). For this purpose, consider the homotopy /i(7,x) from ?/'^(9T) to 
V'|(5T),i.e., 

/i(7,x) = 7V'#(x) + (1 - 7)V'#(x) for 7 G [0, 1]. 

We define 

P = 4(r)-/i#([0,l] xV^^(9T)). (19) 

P is a d-current whose boundary dP is contained in the {d — l)-skeleton of K. Define Q2 = /i#([0, 1] x 
i/j^{dT)). Using the homotopy formula, we get 

dP = d{^liT)-h#{[0,l]x^lidT))) 

= 4(9r)-9/i#([o,i] xv^^(9r)) 

= V#(5T) C (d - l)-skeleton of K. 
Equation (O gives ip)^{T) = P + Q2. Defining Q = Qi + Q2, Equation ^ gives 

T-{P + Q2) = dR + Qi, hence 
T-P = dR + Q. 

Finally, we apply the bounds on the retraction described in Lemma lS!9] and the paragraph following this 
Lemma to the masses of the pushforwards. Noticing that D((t) < A for all a £ K, we get the following 
bounds, which finish the proof of our simplicial deformation theorem (Theorem 15. It . 

M(P) < (4^A-)^' M(r) + A(4^A-)'=+^ M{dT) 

= {4dKf (M(T) + A{AdK) M{dT)) , 
M{dP) < (AdK^'^^MidT), 
M{R) < AM(V'#(r)) 

< A(4^i^)''M(r), and 
M(Q) < AiAdx)'" M{dT) + A(4^a')''"*"^ M{dT) 

□ 

Remark 5. 10. The influence of Simplicial regularity as measured by ki and K2 is clearly revealed by the 
statement of our deformation theorem (Theorem 15. 11 1. Explicit constants are a simple yet useful part of the 
result; as observed above in Remark 15.21 the statement of this theorem leads to an easy observation that the 
flat norm distance between T and P can be made a small as desired by subdividing the simplicial complex 
without increasing the bound on regularity. 

Remark 5.11. We did not explicitly discuss the case of 0-dimensional currents. In this case, the bounds on 
mass expansion are all equal to one. 
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5.2 Comparison of Bounds of Approximation 

Sullivan studied the deformation of integral currents on to the skeleton of a cell complex, which is composed 
of compact convex sets. He presented a deformation theorem for deforming integral currents on to the 
boundary of a cell complex (22. Theorem 4.5]. For ease of comparison, we use our notation to restate the 
bounds given by Sullivan for deforming a d-current T to a polyhedral current P in the boundary of a cell 
complex in W^. Recall that in our simplicial deformation theorem (Theorem 15.11 ). the simplicial complex 
considered has dimension p, and is embedded in for q > p. 




M{P) < rj 2d — — K2 M(r), (20) 



-d+l 

M{dP) < ( )l2tf(^^^K2) ) M{dT), and (21) 




¥{T,P) < {q-d+l)A{M{P)+M{dP)). (22) 

Our results corresponding to the first two bounds in Equations (l20l ) and (|2TI) are presented in Equations (ITOl ) 
and ([TT]) in Theorem 15.11 We use k = p — d to denote the codimension, and thus k = p — d<q — d. 
Hence the exponents of K2 terms in these bounds of Sullivan are at least {k + l)d, while those in our bounds 
are only A; or A; + 1. To obtain the flat norm distance corresponding to the third bound given by Sullivan in 
Equation ((22l) . we use the definition of flat norm distance between two currents specified in Equation (|2]i. 
Using T — P = dQ + R, we combine two of our bounds specified in Equations (fT2l ) and (fT3l) to get 

¥{T,P) < A{4dK)'' {M{T) + {l + 4dK)M{dT)). 

Once again, the exponents of regularity terms ki and K2 in our bound are at most k + 1, while those in 
SuUivan's bound are at least {k + l)d. 

The differences in bounds highlighted above become stark when we are studying high-dimensional 
currents with high codimension. We get much smaller bounds because we deform in a different way. Also, 
we work with simplicial complexes while Sullivan uses the more general cell complexes. At the same time, 
efficient algorithms are more naturally constructed for simplicial complexes than for cell complexes. 



6 Computational Results 

We illustrate computations of MSFN by describing the flat norm decompositions of a 2-manifold with 
boundary embedded in (see Figure 16]). The input set has the underlying shape of a pyramid, to which 
several peaks and troughs of varying scale, as well as random noise, have been added. We model this set 
as a piecewise linear 2-manifold with boundary, and find a triangulation of the same as a subcomplex of a 
tetrahedralization of the 2x2x2 cube centered at the origin, within which the set is located. We use the 
method of constrained Delaunay tetrahedralization 11211 implemented in the package TetGen |[T9]| for this 
purpose. We then compute the MSFN decomposition of the input set at various scale (A) values. At high 
values, e.g., when A = 6, the optimal decomposition resembles the input set with the small kinks due to 
random noise smoothed out. At the other end, for A = 0.01, the optimal decomposition resembles a flat 
"sheet". For intermediate values of A, the optimal decomposition captures features of the input set at varying 
scales. 
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Figure 6: Top left: A view of original pyramidal surface in three dimensions. The remaining three figures 
show the flat norm decomposition for scales A = 6 (top right), A = 2 (bottom left), and A = 0.01 (bottom 
right). See text for further explanation. The images were generated using the package TetView ll20l . 
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The entire 3-complex mesh modehng the cube in question consisted of 14,002 tetrahedra and 28,844 
triangles. For each A, computation of MSFN described above took only a few minutes on a regular PC using 
standard functions from MATLAB. This example demonstrates the feasibility of efficiently computing flat 
norm decompositions of large datasets in high dimensions, for the purposes of denoising or to recover scale 
information of the data. 

7 Discussion 

Our result on simplicial deformation (Theorem I5.lt places the definition of MSFN into clear context. If a 
current lives in the underlying space of a simplicial complex, we can deform it to be a simplicial current 
on the simplicial complex, and do so with controlled error. In fact, by subdividing the simplicial complex 
carefully, we can move this error as close to zero as we like. Since MSFN could be computed efficiently 
when the simplicial complex does not have relative torsion, one could naturally use our approach to compute 
the flat norm of a large majority of currents in arbitrarily large dimensions. 

The MSFN problem, similar to the recent results on the OBCP [8 ], apply notions from algebraic topology 
and discrete optimization to problems from geometric measure theory such as flat norm of currents and 
area-minimizing hypersurfaces. What other classes of problems from the broader area of geometric analysis 
could we tackle using similar approaches? One such question appears to be the following. Is the flat norm 
decomposition of an integral current always another integral current? If not, could we specify conditions 
under which this result holds? Working in the setting of simplicial complexes, results on the existence of 
integral optimal solutions for instances of ILPs with integer right-hand side vectors may prove useful in 
answering this question. 

While L^TV and flat norm computations have been used widely on data in two dimensions, such as 
images, the MSFN opens up the possibility of utilizing flat norm computations for higher dimensional data. 
Similar to the flat norm-based signatures for distinguishing shapes in two dimensions [25], could we define 
shape signatures using MSFN computations to characterize the geometry of sets in arbitrary dimensions? 
The sequence of optimal MSFN decompositions of a given set for varying values of the scale parameter 
A captures all the scale information of its geometry. Could we represent all this information in a compact 
manner, for instance, in the form of a barcode? 
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